Density functional theory and Molecular Dynamics Studies on Energetics and 
Kinetics for Electro- Active Polymers: PVDF and P(VDF-TrFE) 
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We use first principles methods to study static and dynamical mechanical properties of the ferro- 
electric polymer Poly(vinylidene fluoride) (PVDF) and its copolymer with trifluoro ethylene (TrFE). 
We use density functional theory [within the generalized gradient approximation (DFT-GGA)] to 
calculate structures and energetics for various crystalline phases for PVDF and P(VDF-TrFE). We 
find that the lowest energy phase for PVDF is a non-polar crystal with a combination of trans (T) 
and gauche (G) bonds; in the case of the copolymer the role of the extra (bulkier) F atoms is to 
stabilize T bonds. This leads to the higher crystallinity and piezoelectricity observed experimen- 
tally. Using the MSXX first principles-based force field (FF) with molecular dynamics (MD), we 
find that the energy barrier necessary to nucleate a kink (gauche pairs separated by trans bonds) 
in an all-T crystal is much lower (14.9 kcal/mol) in P(VDF-TrFE) copolymer than in PVDF (24.8 
kcal/mol). This correlates with the observation that the polar phase of the copolymer exhibits a 
solid-solid a transition to a non-polar phase under heating while PVDF directly melts. We also 
studied the mobility of an interface between a polar and non-polar phases under uniaxial stress; 
we find a lower threshold stress and a higher mobility in the copolymer as compared with PVDF. 
Finally, considering plastic deformation under applied shear, we find that the chains for P(VDF- 
TrFE) have a very low resistance to sliding, particularly along the chain direction. The atomistic 
characterization of these "unit mechanisms" provides essential input to mesoscopic or macroscopic 
models of electro-active polymers. 
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I. INTRODUCTION 

Poly(vinylidene fluoride) (PVDF) and its copolymers 
with trifluoro ethylene (TrFE) exhibit excellent electro- 
mechanical properties such as ferroelectricity, piezoelec- 
tricity, pyroelectricity, and nonlinear optical properties 
0. Ferroelectricity in PVDF and P(VDF-TrFE) ran- 
dom copolymers was demonstrated with various exper- 
imental techniques [including X-Rays, IR, polarization- 
field hysteresis loops, and the existence of a Curie point 
(only in copolymers)] in the late 70's and early 80's 0. 
Very recently, Zhang and collaborators showed that 
it is possible to make use of electric field-induced phase 
transformations between polar and non-polar phases in 
nano-structured (via electron irradiation) P(VDF-TrFE) 
to obtain large electrostrictive strains (~ 5 %) at high 
frequencies (1000 KHz) and with good energy densities 
(comparable to the best piezoceramics) . Such materials 
are very attractive for many applications requiring soft 
transducers, due to good acoustic impedance match with 
biological tissue and water. 

Molecular dynamics (MD) has previously been used 
to study electro-mechanical properties of PVDF 0, Q 
and P(VDF-TrFE) ^, I, ^ 1, iO , ii J including the pre- 
diction of elastic, dielectric, and piezoelectric constants, 
and polar to non-polar phase transition. However, de- 
spite experimental and theoretical efforts and the tech- 
nological importance of this class of materials, the funda- 
mental molecular processes responsible for their macro- 



scopic electromechanical properties are largely unknown. 
A quantitative characterization of such processes is nec- 
essary to build predictive, first principles-based materials 
models and should be helpful for the design of materials 
with improved properties (such as actuators with higher 
strain, precision, frequency, and energy density). 

In this paper we use atomistic modeling [DFT and first 
principles-based force fields (FFs) with molecular dynam- 
ics (MD)] (described in section II) to study: i) structures 
and energies of various crystal structures (section HI) ; ii) 
the nucleation energy associated with the polar to non- 
polar phase transformation and phase boundary propa- 
gation as a function of applied external stress (section 
IV), and iii) chain sliding under shear deformation (sec- 
tion V). Finally, conclusions are drawn in section VI. The 
comparison between the behavior of PVDF and P ( VDF- 
TrFE) random copolymers allow us to quantify the effect 
of the trifluoro ethylene (TrFE) segments on various ma- 
terials properties. 



A. Crystal structures of PVDF 

Four crystalline polymorphs of PVDF are well charac- 
terized experimentally; they are generally referred to as 
I, II, III, and IV. Following ^ we refer to these as fol- 
lows: I = Tp, 11= TGad, III - TgGpu, and IV = TGpd 
[see Fig. ([Pi]. In this notation we indicate the chain con- 
formation with capital letters (T means all-T, TG means 
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TGTG' and T3G indicates TTTGTTTG') the subscripts 
p or a indicate polar phases with parallel dipoles and 
non-polar phases with anti-parallel dipole moments re- 
spectively; finally the subscripts u and d indicate up- 
up or up-down relative directions of adjacent chains 
When PVDF is cooled from the melt, the TGad phase is 
obtained. The polymer chains contain T and G bonds 
in a TGTG' sequence. This phase has a polar coun- 
terpart: the TGpd phase. The main difference between 
them is the orientation of the dipole moments. In TGpd 
they are parallel while in TGad they are anti-parallel. 
The TGad phase can be converted into TGpd simply by 
poling with an electric field of ~ 100 MV/m. If TGad 
phase is annealed at high-temperature, the polar TsGpu 
phase is formed with a TTTGTTTG' conformation. In 
analogy to the TG phases, there might exist an non- 
polar TsGau phase as sug^sted by Lovinger and 
confirmed by simulations 4r|. In fact, we will use this 
phase to model the non-polar phase with T3G confor- 
mation observed experimentally in electron irradiated 
P(VDF-TrFE). The most interesting polymorph for fer- 
roelectricity is phase Tp; it can be obtained by mechani- 
cal drawing from TsGpu phase or poling from TGpd with 
a ~ 500 MV/m electric field. The polymer chains are in 
an all-T configuration and packed in a parallel fashion. 

Changes in molecular conformation (T, TG, T3G) lead 
to significant shape changes. As shown in Fig. the all- 
T crystal is longer along the chain direction and shorter 
in the directions perpendicular to it than TG and T3G 
phases. This type of structural change is quite differ- 
ent from that in the ceramic ferroelectric compounds in 
which small ionic groups change their dipole orientation 
by rotation and/or displacement. In the polymer sys- 
tem the dipoles are linked together by strong covalent 
bonds and so that orientational change in the dipole mo- 
ments requires cooperative motion of neighboring groups 
through large-scale T-G conformational changes. 



II. METHODS 

The ab initio QM calculations in this paper were per- 
formed using the DFT pseudopotential code SeqQuest 
[T3L IT^ which uses Gaussian basis sets. For all calcu- 
lations (described in Section III) we used the Perdew, 
Burke, and Ernzerhof implementation of the General- 
ized Gradient Approximation (GGA) The SeqQuest 
code calculates atomic forces and the stress tensor which 
were used to relax positions and cell parameters. 

The characterization of many important materials 
properties require the simulation of large systems (thou- 
sands of atoms) for relatively long times (nanoseconds) 
making QM methods impractical. Thus we use the 
First Principles based force field MSXX 0] with MD 
to study energetics, nucleation energies, interface mo- 
bility, and viscoelastic properties. The MSXX force 
field describes the atomic interactions with three energy 
terms: i) electrostatic interactions using QM-derived 



charges; ii) covalent interactions (bonds, angles, tor- 
sions, and cross terms obtained using the Hessian biased 
method to fit QM vibrational frequencies; and iii) van 
der Waals interactions parameterized to reproduce me- 
chanical properties of polyethylene, graphite, CF4 and 
poly(tetrafluoroethylene) crystals. In order to describe 
the TrFE segments we extended the MSXX FF by de- 
scribing the fluorine atom in the CHF group as the F 
atoms in VDF (same atom type and charge), the corre- 
sponding carbon is treated as the CF2 carbon in MSXX 
[denoted Csvf as in with its charge modified to main- 
tain charge neutrality. New three- and four-body terms 
are calculated using combination rules. The supplemen- 
tary material contains all the force field parameters . 
This simple extension of the MSXX FF allows us to char- 
acterize the effects of the presence of TrFE in PVDF. 

All the simulations in this work use periodic bound- 
ary conditions, with the z axis is oriented parallel to the 
polymer chain direction and the y axis in the direction 
of the polarization. 

III. ENERGETICS OF VARIOUS CRYSTALLINE 
PHASES: DFT-GGA AND MSXX FORCE FIELD 

A key issue to understand the electro-mechanical prop- 
erties of these polymers is the relative energies and struc- 
tural changes among the various phases and the origins 
of such differences. We used DFT-GGA and the MSXX 
force field to study the relative energetics and structures 
of various crystalline phases (T, TGad, TGpd, T3Gpu, 
T3Gau) for PVDF and P(VDF-TrFE) (50-50) copolymer 
(50 mole percent VDF). Table © summarizes the ener- 
getics calculated by optimizing atomic positions and cell 
parameters of various phases and isolated chains with pe- 
riodic boundary conditions. Table (|lj also gives lengths 
corresponding to 4 monomers of infinite (periodic) iso- 
lated single chains (in practical terms we performed these 
calculations using three dimensional periodic cells with 
large (50 A) lattice parameters in the directions perpen- 
dicular to the chain direction and optimized the lattice 
parameter along the chain). Here, the crystal cohesive 
energy (Ecoh) is defined as: 

Ecoh = Ecrystal — N X Eisolatod-chainj (1) 

where N is the number of chains inside the crystal 
phase, Ecrystal is the energy of the crystal phase and 
Eisoiatcd-'chain is the energy of corresponding (same chain 
conformation) isolated chain. All energies in Table (P) 
are given in kcal/mol per carbon atom. DFT-GGA pre- 
dicts the TG conformation as the ground state for the 
isolated chains. Next, in relative stability (0.97 kcal/mol 
higher) is the T3G conformation while the all-T chain has 
a much higher energy (1.41 kcal/mol higher than T3G). 
We find that adding fluorine to make P( VDF- TrFE) leads 
to substantial stabilization of the all-T configuration. Al- 
though the TG conformation is still most stable (by 1.08 
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kcal/mol), the T3G structure is higher in energy than 
all-T by 0.06 kcal/mol. Thus the extra F in the copoly- 
mer favors all-T conformations. The MSXX force field 
correctly describes the relative stability among different 
conformations in PVDF and the role of F in P(VDF- 
TrFE) (50-50), in addition to a good description of the 
lengths along the chain direction for various conforma- 
tions. 

In the case of condensed-phase PVDF, DFT-GGA 
yields a large stabilization of the all-T crystal structure 
with a cohesive energy of 2.10 kcal/mol per carbon, fol- 
lowed by T3G phases (Ecoh ~ 1-1 kcal/mol). Finally TG 
phases are less favored by crystallization with Ecoh ^ 0.3 
kcal/mol. This stabilization of T conformations due to 
better packing, represented by the difference in cohesive 
energies, is not enough to reverse the relative energetic 
stability of the phases: our calculations indicate that 
TGpd phase is the most stable crystal followed by the 
TGad phase. The MSXX force field overestimates the co- 
hesive energies, but exhibit trends from PVDF to copoly- 
mer in agreement with QM. This overestimation of Ecoh 
leads to an incorrect ordering of the phases: MSXX pre- 
dicts the all-T phase to be the most stable one. Adding 
F to make P(VDF-TrFE) (50-50) has the same effect as 
in the chains: it favors all-T conformations. QM pre- 
dicts the TsGau phase to be the ground state, followed by 
TsGpu) (0.13 kcal/mol higher) and all-T (0.30 kcal/mol 
higher than TaGau))- Again the effect of F is to favor 
T configurations; this could be an important factor in 
the better crystallinity (and consequently better piezo- 
electric properties) in the copolymer as compared with 
PVDF observed experimentally . 

Zhang et al Q observed experimentally a phase change 
from a polar all-T to non-polar T3G phase when P(VDF- 
TrFE) is irradiated with high-energy electrons. This 
treatment leads to a lower degree of crystallinity and 
consequently shorter coherence length of the crystalline 
regions; this may lead to an increase in the chain-chain 
separation distance. Our simulations explored two ex- 
treme cases of chain-chain separations: perfect crystals 
and isolated chains; our results are consistent with the 
observation in Zhang's experiments and quantify the ef- 
fect of inter-chain cohesive energy. 



IV. PHASE TRANSFORMATION: 
NUCLEATION AND PROPAGATION 

Phase transitions can be induced in PVDF and its 
copolymer with TrFE by temperature, stress or an ex- 
ternal electric field. Understanding the molecular level 
mechanisms responsible for these phase changes is im- 
portant for manufacturing and processing 0] as well as 
to understand the excellent electromechanical properties 
of electron irradiated P(VDF-TrFE) samples jlSj] . 

X-ray diffraction, Fourier transform infrared spec- 
troscopy (FTIR), and differential scanning calorime- 
try have been used to characterize the micro-structural 



changes induced by high energy electron irradiation in 
P(VDF-TrFE) ^18j|. It was found that irradiation reduces 
the degree of crystallinity of the sample decreasing the 
size of the polar domains below a critical value where the 
crystalline regions transform to a non-polar combination 
of T and G bonds (mainly formed by T3G segments). 
This nano-structural change leads to a desirable decrease 
in polarization hysteresis; however for irradiation doses 
higher than the optimum (^^75 Mrad) the hysteresis in- 
creases again |0. This reentrant hysteresis is believed 
to be caused by a high density of cross-linking due to 
high irradiation doses |l8j| . The complex nano-structure 
resulting from irradiation shows excellent electromechan- 
ical properties such as ~5% electrostrictive strain result- 
ing from an electric field induced, reversible phase transi- 
tion between a non-polar T3G and the polar all-T crystal 
0, ITM IT9I I . An atomistic simulations that fully captures 
the nano-structure of the irradiated copolymers would in- 
volve system sizes well beyond current capabilities. Thus 
in this section we characterize two "unit" mechanisms 
that play a key role in the properties of the irradiated 
copolymers and shed some light into phase transitions in 
PVDF and P(VDF-TrFE) in general. 

Previous theoretical studies have focused on phase 
transitions in PVDF and its copolymers, including the 
thermal induced ferroelectric transition % 13, Tlf and 
field induced TGad-TGpd transitions j2Q|, (21,j. These 
phase transitions are generally believed to occur via a 
nucleation and growth process [^l, governed by coopera- 
tive chain rotations and the motion of phase boundaries 

Here we characterize : i) the nucleation of a kink (G 
bond pairs separated by T bonds) in a perfect Tp crys- 
tal; and ii) the mobility of an interface between TsGad 
and Tp phases [the phase transition believed to be re- 
sponsible for the large electrostrictive strain in electron 
irradiated P(VDF-TrFE) To]]. We focus on the differ- 
ence in behavior between PVDF and its 50-50 copolymer 
with TrFE. 



A. Nucleation of a kink in an all-T configuration 

In phase transitions originating from the all-T ferro- 
electric (Ip) phase [such as the ferroelectric-paraelectric 
(F-P PT)], the first unit process is the formation of 
gauche pairs separated by one or more trans bond (de- 
noted as kinks). We use the MSXX FF to calculate the 
energy barrier and molecular level mechanisms associated 
with this process. To do this we start from a perfect Tp 
crystal and impose a harmonic restraint to one dihedral 
angle; we then minimize the total energy of the system 
(including the restraint) as the equilibrium angle of the 
restraint is changed from T to G in small steps (5°). For 
each restraint angle (6*0) we fully relax the system includ- 
ing lattice parameters at zero applied stress. The PVDF 
simulation cell used for this calculation was built from 
the Tp unit cell (two monomers) replicating it 4 times 
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in the x direction, 4 times in y and 20 times in z. The 
resulting cell contains 16 infinite chains (40 monomers 
in each simulation cell per chain) and 1920 atoms. The 
copolymer was obtained from the PVDF cell by randomly 
converting 50 % of the CH2 groups into CHF. 

Fig. ^ shows the total energy as a function of as 
the chosen torsion bond changes from T to G and back to 
T both for PVDF and P(VDF-TrFE) (50-50). For PVDF 
we obtain an energy barrier of 24.8 kcal/mol while for 
the copolymer the barrier reduces to 14.9 kcal/mol. Thus 
our simulations indicate that nucleation of the F-P PT 
takes place much easier in P(VDF-TrFE) than in PVDF. 
Experiments show that under heating PVDF in the Tp 
phase melts without a F-P PT, but in P(VDF-TrFE) the 
F-P PT is observed before melting for concentrations of 
TrFE larger than 18 %. While the molecular level mech- 
anism of the phase transition is not fully understood and 
defects [for instance the interface between crystalline and 
amorphous regions ,21i, ,2^] can play a key role, our cal- 
culated nucleation barriers and the experimental results 
on phase transitions support nucleation and growth as 
the process responsible of the phase transition. 

Since the polymer chains are infinitely periodic (each 
chain is bonded to itself at the cell boundary) a T-to- 
G change in any one dihedral angle must be accompa- 
nied by compensating changes in the opposite direction 
in other bonds. To investigate this effect, we examined 
the change in dihedral angles for bonds neighboring the 
restrained one (denoted as 0) as a function of 6*0, see Fig. 
©. We find that the angles centered around bonds that 
are first nearest neighbors to the restrained one (denoted 
±1) remain T, no change during the process. This is 
because neighboring G bonds are quite unfavorable ener- 
getically. We find instead that it is torsions around the 
second nearest neighbor to the restraint bonds that com- 
pletely compensates for the change in dihedral angle of 
bond 0. Thus as bond transforms to G one of the sec- 
ond nearest neighbors goes to G. The average between 
angles 0, -2, and 2 remains almost constant at the T value 
during the transformation (see magenta line in Fig. |2J|. 

Finally, there is a correlation between the sudden en- 
ergy drop seen in Fig. [3 and abrupt changes in the +2 
and -2 dihedral angles. For small deformations starting 
from the all-T conformation the +2 and -2 angles counter 
the restraint in equal amounts (see Fig. For 9o—90 
° the angle -2 drops back to ^ 180 ° and +2 doubles its 
contribution. This structural relaxation leads to a drop 
in the energy (see Fig. 12)). A similar phenomena occurs 
when the restraint angle is changed back to 180 °. 



B. Mobility of the interface between polar and 
non-polar regions 

We now turn to the second component in the nucle- 
ation and growth process: the velocity at which the al- 
ready nucleated phase transition propagates; this is a 
critical property key to the performance of electron irra- 



diated P(VDF-TrFE). 

Fourier transform infrared spectroscopy shows that 
electron irradiation of P(VDF-TrFE) random copolymers 
decreases all-T content and increases T3G segments T^; 
an electric field-induced phase transition between these 
two phases is responsible for its electromechanical prop- 
erties. Therefore, we model the irradiated material by 
a simulation cell containing non-polar TaGau and all-T 
phases and studied the mobility of the interface between 
them under tensile stress. The interface plane is perpen- 
dicular to the direction of the infinite chains [see Fig. . 
The simulations cell contain 2112 atoms, 4 by 4 chains 
in the x and y directions with each chain containing 44 
monomers in the periodic cell. We equilibrate the sys- 
tem using NPT MD simulations with the MSXX FF at 
T=300 K and zero stress. Using average lattice parame- 
ters obtained from the diagonal components of the shape 
matrix during the NPT run and 90 ° cell angles, we per- 
formed 40 ps of NVT simulations to further equilibrate 
our system. After this procedure residual stress was less 
than 0.05 GPa. The interface between all-T and T3G 
phases did not move during the equilibration process (it 
is well known experimentally that drawing is necessary 
to obtain an all-T configuration). We then applied a ten- 
sile uniaxial stress along the chain direction using NPT 
dynamics; the stress was increased from zero at a load- 
ing rate of 0.5 GPa per ps until the desired stress was 
achieved, then the applied tensile stress was maintained 
and the interface mobility studied. Fig. 10} shows snap- 
shots of this process for P(VDF-TrFE) (50-50) with an 
applied stress of 1.6 GPa. 

Figure [S] shows the percentage of T bonds in the sys- 
tem as a function of time for various applied stresses for 
PVDFlSIa) and P(VDF-TrFE) (50-50) E| (b). Bonds are 
classified as T or G based on the their dihedral angles. 
From the number of G bonds as a function of time we 
compute interface mobility, v, using the following equa- 
tion: 

X L 

where 5n is the change in the number of T bonds during 
the time period 6t normalized to the number chains in 
the simulation cell, and L is the length of the T3G unit 
along chain direction. The calculated interface mobilities 
are plotted in Fig. First, we notice that the thresh- 
old tensile stress to move the interface for P(VDF-TrFE) 
is about 1 GPa, while for PVDF it is much larger {r-^ 2 
GPa). This explains why it is easier to draw non-polar 
phases into polar phases in the case of the copolymer 
as compared with PVDF. Second, we find that the do- 
main interface mobility is much larger in the case of the 
copolymer; for example, we find a velocity of 70 m/s for 
PVDF under a tensile stress of 3 GPa and 100 m/s for 
P(VDF-TrFE) under 2.8 GPa. These velocities are less 
than onc-hundrcdth of the speed of sound propagating 
in PVDF and P(VDF-TrFE) [The sound velocity in the 
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<001> direction is about 12 km/s and 11 km/s for PVDF 
and P(VDF-TrFE) Tp phases respectively.] We note that 
MD simulations incorporate in a natural way the many- 
body effects and cooperative phenomena involved in the 
motion of a domain boundary. 

V. CHAIN SLIDING UNDER SHEAR 
DEFORMATION 

The various crystalline structures of PVDF and 
P(VDF-TrFE) are formed by covalently bonded chains 
packed in a parallel fashion interacting via weak van der 
Waals forces. Since the various phases have large differ- 
ences in length along the chain direction, we expect inter- 
chain slip to be an important mechanisms during the 
phase transition, especially when the interface separat- 
ing the two phases is along the chains. This process may 
play an important role to release deviatoric stresses that 
build up during phase transitions [as when electrostric- 
tive P(VDF-TrFE) is used as an actuator]. 

Two popular methods for studying the viscous prop- 
erties in molecular simulations are: i) the Green-Kubo 
(GK) method to analyze the time integral of the stress- 
stress correlation function during equilibrium molecular 
dynamics leadin g d irectly to the ordinary viscosity at 
zero strain rate |23l ]. and ii) non-equilibrium molecular 
dynamics (NEMD) in which a deformation is imposed on 
the system and the time-averaged resulting component of 
the stress tensor is calculated; the stress-strain rate data 
can then be used to calculate viscosity as a function of 
strain rate. In this paper we use the NEMD technique, 
imposing volume-conserving pure shear deformation to 
our simulation cells in two directions for various strain 
rates. 

Using the convention that the c axis is the chain di- 
rection and the b axis is parallel to the dipole moments 
of the Tp phase, we deformed Tp P(VDF-TrFE) crystals 
along the following systems: i) on the (100) plane in the 
(010) direction (denoted as perpendicular to the chains), 
and ii) on the (010) plane in the (001) direction (denoted 
as along the chains). We applied strain rates in the range 
1.66 X 10^° 1/s to 6.67 x 10^*^ 1/s using simulations cells 
containing 2x4x4 unit cells (384 atoms). Symbols 
in Fig. O show the resulting averaged shear stress as 
a function of the imposed strain rate for both the per- 
pendicular and parallel cases. We fitted our data to the 
widely used inverse hyperbolic sine flow equation: 

.(^) = ^|Isinh- (1) , (3) 

where a is the shear stress, 7 is the strain rate, Vo de- 
notes the effective volume of the unit deformation, and 
7o is a characteristic frequency (inverse time). The fit- 
ted functions are shown as solid lines in Fig. (0. The 
fiow equation [Eq. Q] is very general in nature and can 
describe various mechanisms such as free volume-based 



deformation in metallic glasse s l24l l25l and dislocation- 
based single crystal plasticity |26l l27j . The resulting ef- 
fective volume in the parallel case is Vq = 281. 4A'^ and 
in the perpendicular direction is 116.3 A^] the charac- 
teristic frequencies 70 are: 1.33 x 10^°l/s and 0.17 x 
10^" 1/s in the parallel and perpendicular directions re- 
spectively. As expected, we find that the direction of 
easy slip has a larger characteristic frequency [70 is pro- 
portional to exp (AG/kfiT) where AG is an activation 
barrier p^l25ll26| ] and a larger effective volume. 

We calculate viscosity (77) as 77 = cr/7; Figure ||SJ shows 
the viscosity as a function of strain rate for both direc- 
tions. The symbols represent the atomistic data and the 
lines were obtained from the previous fits. We see from 
Fig. (jHJ that we are in the shear-thinning regime (an 
increase in the shear rate leads to a decrease in viscos- 
ity) and the viscosity along the chains direction (the one 
important to release stress during a phase transition) is 
significantly smaller than that in the perpendicular di- 
rection. This is somewhat evident by examining the Tp 
crystal structure and more clear by analyzing the MD 
trajectories. Fig. (9) shows snapshots from our runs for 
7 = 4.0 x 10^" both in the parallel (top panels) and per- 
pendicular (bottom panels) directions. While both defor- 
mations lead to disorder of the chains, the perpendicular 
case involves much larger molecular rotations and rela- 
tive shifts, and consequently larger energies and stresses. 



VI. CONCLUSIONS 

We used a combination of ab initio (DFT-GGA) cal- 
culations and MD with First Principles Force Fields to 
characterize several "unit mechanisms" that govern the 
electromechanical properties of elcctroactive polymers 
PVDF and P(VDF-TrFE) (50-50). The quantum me- 
chanics calculations of energetics and structures of var- 
ious phases of PVDF and P(VDF-TrFE) show that T 
configurations are stabilized energetically by the addition 
of bulkier F atoms in TrFE. This result is likely related 
with the better crystallinity and piezoelectric properties 
of P(VDF-TrFE), because of its stabihzation of the all-T 
form, Tp. Our results also show that conformations with 
T and G bonds are energetically favorable for large inter- 
chain separations as compared with all-T structures; this 
observation is consistent with the experimental observa- 
tion that samples bombarded with high-energy electrons 
favor T3G non-polar conformations. 

Our MD results show that the energy barrier required 
to nucleate a G bond in an all-T configuration in P(VDF- 
TrFE) is 40% lower than in PVDF. Experimentally 
copolymers show a P-N PT upon heating before melting 
while PVDF directly melts; this is consistent with the 
calculated nucleation barriers. 

We also studied the motion of the interface separat- 
ing non-polar TsG au and polar Tp phases under uniaxial 
stress to model the reversible phase transition responsible 
for electrostriction in electron irradiated P(VDF-TrFE). 
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We found a smaller threshold stress and higher domain 
waU mobility in P(VDF-TrFE) as compared to PVDF. 
Interface mobility is an important property for actuators 
and MD provides a detailed understanding and quantita- 
tive characterization of this process that should be useful 
in mcso- or macro-scopic simulations. Finally, our sim- 
ulations of the dynamics of chain sliding reveal that slip 
occurs more easily along the direction of the chains than 
perpendicular to them. 

In summary we have used First Principles-based multi- 
scale modeling to characterize dynamical and static prop- 
erties PVDF and its copolymers with TrFE. Our simula- 
tions are consistent with available experimental observa- 
tions. In order to make a more quantitative comparison 
between our dynamical calculations and experiments, we 
are currently using the atomistic results with mesoscopic 
modeling ^2^] to predict the behavior of a real actuator 



whose performance can be directly compared to experi- 
ments. We foresee that such multi-scale simulations will 
be an important tool to guide the design of new materials 
with improved properties. 
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TABLE I: Energetics of crystalline phases and isolated infinite chains of PVDF and P(VDF-TrFE) (50-50). All energies are 
given in kcal/mol per carbon atom. The reference energies for crystal phases and chains are the energy of Tp phase and all-T 
chain, respectively. The cohesive energy defined by Eq. (1) in the text is denoted by Ecoh. The length (L) of isolated chain is 
defined as the length of four monomers along the chain direction. (QM) denotes DFT-GGA data; (EE) denotes MSXX force 
field. 





PVDF 


P(VDF-TrFE) 


Crystal 


AE (QM) 


Ecoh(QM) 


AE (EE) 


Ecoh (FF) 


AE (QM) 


Ecoh(QM) 


AE (FF) 


Ecoh (FF) 


Tp 





-2.10 





-3.98 





-0.86 





-3.02 


TGad 


-0.59 


-0.31 


0.39 


-2.40 


1.86 


2.09 


1.17 


-0.72 


TGpd 


-0.62 


-0.34 


0.35 


-2.45 


2.42 


2.65 


1.26 


-0.62 


TaGpu 


-0.51 


-1.20 


0.19 


-3.14 


-0.17 


-1.10 


0.58 


-2.49 


TsGau 


-0.38 


-1.03 


0.23 


-3.10 


-0.30 


-1.23 


0.56 


-2.50 


Single Chains 


AE (QM) 


AE (EE) 


L (A) (QM) 


L {A) (FF) 


AE (QM) 


AE (FF) 


L (A) (QM) 


L (A) (FF) 


T 








5.22 


5.13 








5.19 


5.21 


TG 


-2.38 


-1.19 


4.66 


4.56 


-1.08 


-1.13 


4.71 


4.68 


TsG 


-1.41 


-0.65 


4.62 


4.54 


0.06 


0.05 


4.72 


4.70 








m 


4 4 i 



TTTT(p) TGTG'(a) TGTG'(p) TsGTsG'Cp) T3GT3G'(a) 



FIG. 1: Crystal structures of various phases for PVDF. The polar phases with parallel dipole moments are denoted by (p), 
while non-polar phases with anti-parallel dipole moments are denoted by (a). The existence of the non-polar TaGad phase was 
conjectured by Lovinger ,12| and proved by Karasawa and Goddard Q. 




Torsional angle (deg) 



FIG. 2: Energy barrier to nucleate a G bond in at Tp perfect crystal both PVDF and P(VDF-TrFE) (50-50). We used supercells 
consisting of 16 chains each consisting of 8 monomers (2112 atoms). 
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FIG. 3: Evolution of various torsion angles in PVDF as a function of the restraint dihedral angle as it is transform from T 
to G. The diagram at the top shows the notation: the restraint is applied to torsion angle for the bond denoted 0. We find 
that the first nearest neighbor bonds to the restraint one (+1 and -1) remain in the T conformation, while the second nearest 
neighbor bonds (+2 and -2) change to compensate for the changes in bond 0. The purple curve (pluses) shows that the average 
of the torsion angles remains at 180 degrees, which is topologically necessary since the chains are infinitely periodic. Fig. 
corresponds to restraint varying from 180 to 50 while Fig. reflects the reverse process. 




time = 4 ps = 8 ps 

FIG. 4: Snapshots of the TsGad-Tp interface for P(VDF-TrFE) (50-50) under 1.6 GPa tensile stress as various times obtained 
from the MD simulation. 




FIG. 5: TaGad-Tp interface mobility. Number of T bonds as a function of time for PVDF (Fig. and P(VDF-TrFE) (50-50) 

(Fig. (Et)) under various uniaxial stress states at T=300 K. 
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FIG. 7: Shear stress and a function of applied strain rate for deformation in two different orientations with respect to the 
polymer chains for P(VDF-TrFE) (50-50) at 300K. Lines represent fits to the MD data using the hyperbolic sine function . 




FIG. 8: Viscosity for two different orientations with respect to the chains of P(VDF-TrFE) (50-50) at 300K. Points represent 
MD data and hnes were obtained from the fits to the stress-strain rate data using Eq. 3. 
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(9a) Shear strain applied parallel to chains 
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(9b) Shear strain applied perpendicular to chains 
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FIG. 9: Snapshots obtained form MD simulation of shear deformation of a P(VDF-TrFE) (50-50) crystal using strain rate of 
4.0 X 10^" 1/s in the direction parallel (a) and perpendicular (b) to the polymer chains. 



